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Abstract 

Background/Aim: Efforts have been made to eliminate wild poliovirus transmission since 1988 when the World 
Health Organization began its global eradication campaign. Since then, the incidence of polio has decreased 
significantly. However, serotype 1 and serotype 3 still circulate endemically in Pakistan and Afghanistan. Both 
countries constitute a single epidemiologic block representing one of the three remaining major global reservoirs 
of poliovirus transmission. In this study we used genetic sequence data to investigate transmission links among 
viruses from diverse locations during 2005-2007. 

Methods: In order to find the origins and routes of wild type 1 poliovirus circulation, polioviruses were isolated 
from faecal samples of Acute Flaccid Paralysis (AFP) patients. We used viral cultures, two intratypic differentiation 
methods PGR, ELISA to characterize as vaccine or wild type 1 and nucleic acid sequencing of entire VPl region of 
poliovirus genome to determine the genetic relatedness. 

Results: One hundred eleven wild type 1 poliovirus isolates were subjected to nucleotide sequencing for genetic 
variation study. Considering the 15% divergence of the sequences from Sabin 1, Phylogenetic analysis by MEGA 
software revealed that active inter and intra country transmission of many genetically distinct strains of wild 
poliovirus type 1 belonged to genotype SOAS which is indigenous in this region. By grouping wild type 1 
polioviruses according to nucleotide sequence homology, three distinct clusters A, B and C were obtained with 
multiple chains of transmission together with some silent circulations represented by orphan lineages. 

Conclusion: Our results emphasize that there was a persistent transmission of wild typel polioviruses in Pakistan 
and Afghanistan during 2005-2007. The epidemiologic information provided by the sequence data can contribute 
to the formulation of better strategies for poliomyelitis control to those critical areas, associated with high risk 
population groups which include migrants, internally displaced people, and refugees. The implication of this study 
is to maintain high quality mass immunization with oral polio vaccine (OPV) in order to interrupt chains of virus 
transmission in both countries to endorse substantial progress in Eastern-Mediterranean region. 
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Background 

Since 1988 the world has come very close to eradicate 
polio through global polio eradication initiative [1]. The 
objective of this initiative was to interrupt wild PVs as 
soon as possible to achieve certification of Global Polio 
Eradication and to strengthen the routine immunization 
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and surveillance as well Since its inception, this initia- 
tive has made remarkable headway worldwide and num- 
ber of cases of poliomyelitis dropped from 35,000 in 
1988 to 1650 in 2008 [2]. However, the Indigenous 
transmission of wild poliovirus (WPV) has never been 
interrupted in Afghanistan, Pakistan, India, and Nigeria 
[3-5]. 

The Polio Eradication Initiative was launched in 
Pakistan in 1994 with remarkable success, the number 
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of confirmed cases declined from 1015 in 1997 to 32 
wild polio cases including 19 WPVl and 13 WPV 3 in 
only 18 of 120 districts of Pakistan [6,7]. In Afghanistan 
polio eradication activities started in 1997, [8] since 
then, significant progress have been made in both coun- 
tries, wild type 2 poliovirus has been knocked out from 
Pakistan since 1997 and in Afghanistan 1999 [7,9]. 

Despite intensified mass immunization activities, the 
two countries failed to interrupt poliovirus transmission. 
Persistent pockets of polio transmission along the bor- 
der between Afghanistan and Pakistan are key epidemio- 
logical challenges due to insecurity and continued 
conflict [10,11]. 

Poliovirus, the causative agent of poliomyelitis, is a 
human enterovirus and member of family Picornaviridae 
[12]. It is composed of single- stranded positive- sense RNA 
genome that is about 7500 nucleotides long [13]. There 
are three serotypes of poliovirus, PVl, PV2, and PV3; each 
with a slightly different outer capsid protein which define 
cellular receptor specificity and virus antigenicity. There- 
fore infection with one serotype does not prevent infection 
with another serotype [14,15]. PVl is the most common 
type encountered in nature, however aU three types are 
extremely infectious and can circulate independently 
[16,17]. Long excretion periods and low population immu- 
nity also support its rapid evolution and spread in humans 
[18-20]. The mutation rate in poliovirus is relatively high 
with a synonymous substitution rate of 1.0 x 10'^ substitu- 
tions/site/year [21,22] which creates a wide variety of 
mutants that are referred as different genetic strains or 
genotypes. Each serotype of wild poliovirus has many dif- 
ferent genotypes that are distributed geographically and 
co-circulate worldwide [20,23,24]. Genetic diversity of wild 
polioviruses is determined by its genome sequencing that 
not only provide genetic relationship among wild Polio- 
viruses as well as it helps to monitor the progress of Polio 
Eradication Program [25,26]. 

In current study 111 wild typel poliovirus isolates col- 
lected during 2005-2007 from Pakistan and Afghanistan 
were sequenced for complete VPl gene to study their 
genetic diversity and to reveal the indigenous genotype 
in this region that would be helpful to track the trans- 
mission patterns in Pakistan and Afghanistan being a 
single epidemiological block. 

Methods 

Virus isolation & identification 

Stool samples of children with AFP (Acute Flaccid 
Paralysis) and from apparently healthy contacts were 
collected as per WHO guidelines within 14 days of 
paralysis onset and were sent to WHO Regional Refer- 
ence Laboratory for Polio Eradication Initiative, Depart- 
ment of Virology, National Institute of Health, 
Islamabad, Pakistan. 



These samples were processed for isolation and identifi- 
cation of poliovirus by standard techniques as described in 
World Health Organization manual [27]. Briefly, 200ul of 
stool suspension extracted with chloroform was inoculated 
in culture tubes of two cell lines; RD cell line (derived 
from rabdomayosarcoma cells of humans) and L20B cell 
line (mouse lymphoma cells expressing the human polio- 
virus receptor) [28]. Positive samples showed a character- 
istic poliovirus cytopathic effect (CPE), i.e. rounded, 
refractive cells were serotyped by microneutralizaton using 
pools of antisera against polioviruses as recommended by 
WHO [27,29]. 

Intratypic differentiation 

To differentiate Sabin (SLl) and Wild type (NSLl) 
Polioviruses, PGR and ELISA were used according to 
WHO recommendations [30,31]. 

Genomic sequencing 

The viral RNA of wild type 1 polioviruses was extracted 
from 140 (il of infected tissue culture fluid using the QI 
Aamp viral RNA kit (QIAGEN GmbH, Hilden) and subse- 
quently stored in aliquots at-70°G. Reverse transcription 
Polymerase chain reaction (RT-PGR) was performed in 
single step of 50 (il reaction volume. RT-PGR was used to 
amplif)^ the entire VPl capsid protein coding region of 906 
nucleotides. The primers used in these reactions listed in 
(Table 1). PGR products were purified using the QIAquick 
PGR purification kit (QIAGEN GmbH, Hilden). The puri- 
fied amplicons were further processed in cycle sequencing 
reaction using ABI Prism Big Dye Terminator cycle 
sequencing Ready Reaction kit (dGTP BigDye®, Applied 
Biosystems) with the overlapping primers (Table 1). The 
cycle sequencing program was 25 cycles of 20 s at 94°G, 
15 s at 42°G and 4 mins at 60°G. After the completion of 
program, the reactions were cooled to 4°G. The purified 
product of sequencing reaction was suspended in deio- 
nized formamide and run in an automated Genetic Analy- 
zer (ABI Prism, model 3100). 

Sequence analysis softwares 

Sequence data was analyzed with the help of sequencher 
software (Gene Godes v. 4.5) which provides different 
methods to edit and align sequence fragments into 

Table 1 Sequences of primers used in cDNA synthesis. 



PCR and sequencing in this study 



Primer 


Location (nt) 


Sequence 


Y7 


2396-2418 


5'-GGGmGTGTCAGCCTGTAATGA-3' 


Q8 


3504-3485 


5'-AAGAGGTCTCTRTOCACAT-3 


PVl A * 


2954-2935 


5'-^IAIIGCRTGICCR^R^-3' 


PV4 S * 


2830-2849 


5'-ACITAYAARGAYACIGTICA-3' 



* Inosine Containing Primeers, Use degenerate PCR Condition 
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larger structure (contigs). Phylogenetic analysis of wild 
type 1 polioviruses was carried out by using Molecular 
Evolutionary Genetic Analysis (MEGA v.4) [32]. Evolu- 
tionary distances were calculated using Kimura 2-para- 
meter model [33] and phylogenetic tree was constructed 
using the neighbor-joining method [33,34]. 

Results 

Study isolates 

A total of 25,965 and 4,485 AFP cases including contact 
cases from Pakistan and Afghanistan respectively were 
processed during 3 years time period and 111 wild 
typel polioviruses were detected as shown in (Tables 2 
and 3). 

Genetic comparison of wild type 1 polio isolates 

The sequences of VPl gene for all 111 wild typel polio- 
virus isolates were compared with each other to measure 
their genetic similarity by using Neighbor-Joining with 
Kimura 2 Parameter [33,35]. Based on pairwise distance 
comparison of closely related wild type 1 polioviruses, a 
dendogram of sequence relationship between isolates was 
constructed (Figure 1). In these comparisons, we followed 
imperative factor which defined that viruses with < 1% 
difference from Sabin vaccine virus are classified as Sabin- 
like; those with 1-15% difference as vaccine-derived polio- 
viruses and those with > 15% difference as the wild virus 
[36]. The criteria for polioviruses genotypes has been 
defined as groups of strains that show more than 85% 
sequence homology in the VP 1/2 A region and a cluster is 
defined as a group of isolates showing > 95% sequence 



similarity. A strain represents new cluster if it has > 5% 
nucleotide sequence divergence with previous ones 
[21,37]. 

As a result all 111 isolates were grouped into three 
major clusters A, B, & C respectively. Within these clus- 
ters, a definite pattern of virus distribution was observed 
representing lineages i.e., cluster A with six (a-f) corre- 
sponding lineages (Figure 2) and cluster B isolates have 
three (g-i) different lineages (Figure 3). While the cluster 
C has no further lineage because it has single wild typel 
virus isolate. The true root of tree cannot be explicitly 
determined from the sequence data; however a potential 
mid root Pak 93-5218 was used to connect the sequences 
of all 3 years (2005-2007) isolates, represents a hypotheti- 
cal genetic founder. 

Geographic distribution of wild type 1 polioviruses based 
on genetic data 
Cluster A 

This cluster comprised of 81 wild poliovirus strains. The 
degree of divergence in VPl gene sequence ranged from 
0.6% to 18.7%. In this cluster the virus reservoirs were 
confined among the various geographical areas of both 
countries. Six genetic lineages a, b, c, d, e & f (Figure 2) 
prevail within the cluster A, some remained active over 
the 3 years while the others indicate successful termina- 
tion through effluent control measures. 

Lineage a (Figure 2) was observed with a higher pro- 
portion of wild polioviruses having diversified and vast 
geographical span in various districts of both countries 
(Figure 4). It also reveals that strains from same lineages 



Table 2 Distribution of wild typel polioviruses in Pakistan with total number of districts infected during 2005-2007 









2005 






2006 






2007 


Month 


Wild 
PI 


Province 


No. of infected 
districts 


Wild 
PI 


Province 


No. of infected 
districts 


Wild 
PI 


Province 


No. of infected 
districts 


January 


4 


BN, PB, 
SD 


4 








2 


KP, SD 


2 


February 








2 


KP, BN 


2 








March 


2 


SD, BN 


2 














April 


1 


PB 


1 


1 


SD 


1 








May 


3 


PB 


2 


1 


BN 


1 


1 


KP 


1 


June 


2 


PB, SD 


2 


3 


KP, BN 


3 


2 


KP 


2 


July 


2 


PB, SD 


2 


3 


KP, PB 


2 








August 


3 


KP.BN 


3 


2 


KP, BN 


3 


2 


SD 


2 


September 


2 


BN, SD 


2 


3 


KP, SD 


3 


2 


SD 


1 


October 


2 


PB, KP 


2 


2 


KP 


2 








November 


3 


PB, KP, 
BN 


3 


2 


SD 


2 


3 


KP, BN 


3 


December 


2 


BN, SD 


2 








8 


KP, SD 


4 


Total (AFP 
+Contact) 


26 




25 


19 




19 


20 




15 
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Table 3 Distribution of wild type1 polioviruses in Afghanistan with total number of districts infected during 
2005-2007 







2005 




2006 






2007 


Month 


Wild 
PI 


Province No. of infected 
districts 


Wild 
PI 


Province 


No. of infected 
districts 


Wild 
PI 


Province No. of infected 
districts 


January 






1 


KANDHAR 


1 






February 






3 


KANDHAR 


2 






March 






2 


KANDHAR 


2 






April 






1 


KANDHAR 


1 


2 


HEL+KAND 2 


May 






7 


FARAH+HEL+KAND 
+ORU+ZABOL 


5 


1 


LAGHMAN 1 


June 






9 


FARAH+HEL+KAND 
+ORU 


10 


1 


NANGARHAR 1 


July 






3 


KANDHAR 


4 






August - - - - - _ _ _ 


September 


2 


NAGARHAR 1 


3 


NANG+ORU 


3 






October 


1 


HELMAND 1 






2 






November 


1 


HELMAND 1 


1 


BAGHLAN 


1 






December 


2 


HEL 2 
+KAND 








2 


HELMAND 2 


Total (AFP 
+Contact) 


6 


5 


33 




31 


6 


6 



remained circulating throughout 3 years time period. 
These strains were isolated in different districts during 
2005 from Pakistan such as one district of Sindh; 
[Larkana (SD/51/05/007)] two districts of Balochistan; 
[Quetta (BN/11/05/021)] and [Pishin (BN/12/05/007)] 
one from district of KP [Tank (NW/63/05/006)]. On the 
other hand one of the isolate AFG/08/05/456 reported 
from District Kandhar in Afghanistan indicates the 
south-western circulation during 2005 between both 
countries. The virus descendants from the same lineage 
were introduced into different infected districts of Paki- 
stan as well as Afghanistan during the next 2 years 2006 
and 2007. 

A massive wild type 1 polio virus circulation was also 
observed in Afghanistan during 2006, such as one isolate 
AFG//06/06/494 from Province (Nangarhar) of district 
Rodat, one isolate AFG/01/06/562 from (Baghlan) dis- 
trict Puli Khumri, three isolates AFG/08/06/508,541,543 
from (Helmand) district Bust and two isolates from this 
lineage AFG/08/06/461,532 from kandhar with districts 
SpinBoldak and Maywand. Similarly in Pakistan during 
same period, polio wild type 1 virus circulation was seen 
in Killa Abdullah and Quetta districts (Balochistan), 
Bajour district (KP) and Karachi Korangi district (Sindh) 
while in 2007, Kambar & Thatta (Sindh); Lasbela (Balo- 
chistan); Nowshera, South Waziristan & Khyber (KP) 
were infected. However, no strain of this lineage was 
observed during 2005-2007 in Punjab. 

The lineage 'b' (Figure 2) remained confined to Punjab 
province of Pakistan, infecting districts Dera Ghazi Khan 



and Khanewal in 2005, while in 2006 a single case of 
poliovirus paralysis was identified from district Muzaf- 
fargarh. At the same time, two viruses were also docu- 
mented from Afghanistan (Kandhar and Zabul) that 
again supports an unseen intermediate region of virus 
circulation between two countries. During 2007 this 
lineage remained silent and no virus was detected from 
both countries. 

The circulation of lineage and lineage'd' (Figure 2) 
was observed in a very limited districts of Punjab and 
Sindh and also been found in only Helmand province of 
Afghanistan during 2005-2006 which is closely related 
to the virus isolated from Jaffarabad district of Balochi- 
stan province in Pakistan. 

The lineage 'e' (Figure 2) sustained between both coun- 
tries showing diffused circulation within the Balochistan, 
Helmand and Kandhar with subsequent transmission to 
Sindh during 2006 and 2007. However the descendents 
of this lineage are very important for Afghanistan as 
most of the wild type 1 polio isolates showing wide 
spread circulation in Kandhar, Oruzgan, Helmand and 
Farah during 2006. 

It was also evident that adverse infection has been 
persisted in Killa Abdullah district of Balochistan for the 
three consecutive years 2005-2007 contributing to the 
transmission of lineage 'f (Figure 2) during 2005 and 
later on spread to Dera Ismael Khan district of KP. 

None of the AFP case during the year 2007 matched 
with this lineage giving a positive indication of its suc- 
cessful disruption. 
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y— AFG/08/05;456 Kimdhur CS (3( 
|_ I — BN/14/06/005 ICAbdullah (23/2 



j— AFlOOl/06/562 Baghai 



- Are/08/06/523 Oruzgan (18 



Figure 1 A dendogram summarizing the sequence relationship 
among wild typel polioviruses isolated in Pakistan 2005-2007 
based on VP1 nucleotide sequences. Tree was constructed by 
Neighbor-Joining metliod using Kimura2-parameter model. Isolates 
are represented by epid no, date of onset and province codes. 
Codes were mentioned as PAK: Pakistan, AFG: Afghanistan, BN: 
Balochictan, KP: Khyber Pakhtoonkahw, PB: Punjab, SD: Sindh. 



- SD/55/07/042 (6052)0 

— SD/55/07/042 Kamber CI (11/8/2007) 

- BN/55/07/01 1 Lasbclla (23/1 1/2007) 
5/07/028 Kamber (9/1 1/2007) 

— SD/55/07/034 Kamber (30/9/2007) 

— SD/22/07/019 Thatta (7/8/2007) 
- SD/68/06/027 KHLKoraiigi (30/11/2006) 
1 KHLKorangi (14/9/2006) 
- BN/1 1/06/01 1 Quetta (23/5/2006) 

-AFG/0 1/06/562 Baghan (28/1 1/2006) 
- BN/ 1 2/05/007Pis hiii(6082 )1 2/3/2005 

I NW/33/06/030 Bajour (17/9/2006) 

L- AFG/06/06/494 Nangarhar (29/9/2006) 

I NW/34/07/001 Khyber (30/1/2007) 

I NW65/07/01 2 S.Waziristan (2 1 /1 1 /2007) 

I NW/32/07/024 NoKshera (14/5/2007) 

-BN/11/05/021 Qiictta (13/9/2005) 
I — AFG/08/06/543 Helmand(19/6/2006 

H| AFG/08/06/532 Kandliar (8/6/2006) 

j I AFG/08/06/54 1 Hel mand ( 1 5/6/2006 ) 

I AFG/08/06/508 Helmand (24/7/2006) 

BN/14/06/014 KAbduIlah (2/6/2006) 

BN/14/06/005 KAbdiillah (23/2/2006) 

AFC;/08/06/461 Kandhar (4/2/2006) 

I AFG/08/05/456 Kandhar C5 (30/12/2005) 

■ SD/51/05/007 Larkana (21/3/2005) 

NW/63/05/006 Tank (5/6/2005) 

PB/81/05/012 D.G.Khan (1 1/4/2005) 
-PB/8I/05/019 D.G.Khan (10/5/2005) 



— PB/8 1/05/031 D.G.k 

- PB/62/05/001 Khanew 



an (6/6/2005) 
(6/1/2005) 



PB/81/05/062 D.G.Khan (28/11/2005) 




PB/83/06/043 Muzafarghar (21/7/2006) 
'~i pPB/81/05/021 D.G.Khan (20/5/2005) 
H r- PB/81/05/059 D.G.Khan (17/5/2005) 

L AFG/08/06/465 Kandhar (10/2/2006) 

^ AFG/08/06/524 Zabid (17/5/2006) 

AFG/08/05/439 Helmand (1 1/12/2005) 

SD/53/05/082 Jacohabad (30/12/2005) 

- BN/43/05/029 Jafarabad (27/9/2005) 

PB/62/0S/016 Khanenal (21/5/2005) 
PB/72/05/051 R.Y.khan (21/7/2005) 

;/06/003 KHLG.lqbal(28/2/2006) 
SD/32/06/049 Sanghar (21/1 1/2006) 
AFG/08/06/458 Kandhar (18/1/2006) 

SD/61/07/009 KHLBaldia (5/8/2007) 
SD/42/07/049 Ghotki (27/12/2007) 
AFG/08/07/757 Hel mand (25/12/2007) 
AFG/08/07/641 Helmand(10/4/2007) 
AFG/08/06/565 Kandhar (24/7/2006) 
D/61/07/001 KHLBaldia (1/12/2007) 
AFG/08/07/645 Kandhar (16/4/2007) 

BN/14/07/012 KAbdiillah (1 1/1 1/2007) 
AFG/08/07/759(1660) 
p BN/25/05/002 Ksaifullah (10/8/2005) 
-AFG/08/05/425 Helmand (29/10/2005) 
|AFG/03/06/611 Farah (23/5/2006) 
n-AFG/03/06/611 (694)C2 
J L AFG/08/06/487 Kandhar (19/3/2006) 

I AFG/08/06/560 Kandhar (13/9/2006) 

^ AFG/08/06/556 Kandhar C4 (27/7/2006) 

AFG/08/06/5 1 7 Helmand ( 1 0/5/2006) 

I BN/14/06/016K.Abduliah(4070)16/8/2006 

n AFG/08/06/557 Kandhar (10/7/2006) 

I AFG/08/06/579 Oruzgan (7/9/2006) 

|AFG/08/06/490 Kandhar C2 (8/3/2006) 
I I AFG/08/06/490 Kandhar (8/3/2006) 

Ij AFG/08/06/530 Ourzgan C3 

T AFG/08/06/552 Kandhar (17/6/2006) 

AFG/03/06/61 1 (695)0 

AFG/08/05/427 Helmand (6/1 1/2005) 

I — AFC;/08/06/528 Kandhar (26/5/2006) 
J — AFG/08/06/540 Helmand(12/6/2006) 

I AFG/08/06/545 Kandhar (10/6/2006) 

p AFG/08/06/509 Kandhar (1/5/2006) 
|L AFG/08/06/5 18 Kandhar (4/5/2006) 

] I AFG/08/06/523 Oruzgan (18/5/2006) 

I AFG/08/06/581 Oruzgan (2/9/2006) 

AFG/08/06/463 Kandhar (8/2/2006) 

BN/14/05/015 K-Abdidlah (24/11/2005) 
NW/60/06/002 D.LKhan (23/2/2006) 



Figure 2 Dendogram of Cluster A summarizing six lineages 
(a-f) sequence relationship among wild typel polioviruses 
isolated in Pakistan 2005-2007 based on VP1 nucleotide 
sequences. Tree was constructed by Neighbor-Joining method 
using Kimura2-parameter model. Isolates are represented by epid 
no, date of onset and province codes. Codes were mentioned as 
PAK: Pakistan, AFG: Afghanistan, BN: Balochictan, KP: Khyber 
Pakhtoonkahw, PB: Punjab, SD: Sindh. 
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Cluster B 

This cluster contained three lineages g, h, and i (Figure 3) 
each one with a peculiar confined pattern of circulation 
in Afghanistan and Pakistan. The sequence similarity of 
cluster B in VPl gene ranged from 0.1% to 9.1%. The iso- 
lates under lineage 'g' (Figure 3) has a broad geographical 
circulation among three provinces Punjab, Sindh and 
Balochistan during 2005 with no further onward trans- 
mission during 2006 and 2007. 



The lineage 'h' (Figure 3) remained confined only to 
KP province circulating in two districts Peshawar and 
Bajour during 2005. With the undetected circulation 
during 2006, the same lineage re-appeared in Peshawar 
during 2007 (NW/30/07/058) sharing 2.8% genetic simi- 
larity with the isolate (NW/30/05/098) reported in 2005. 
This isolate also showed maximum VPl homology with 
two Afghanistan isolates AFG/06/07/578 (2.8%) and 
AFG/06/07/568 (2.9%) from district Nangarhar and 



— SD/42/05/017 Ghotki (25/7/2005) 

SD/51/05/037 Larkana (29/9/2005) 

-BN/31/05/005 Sibi (27/3/2005) 



Lineage g 



k; 



PB/61/05/005 Multan (29/1/2005) 
BN/24/05/003 Musakhel C2 (22/7/2005) 

ArG/06/07/578 Nangarhar (18/6/2007) 
NW/30/07/058 Peshawar (24/6/2007) 
AFG/06/07/568 Lagman (26/5/2007) 

NW/30/05/098 Peshawar (15/10/2005) 

NW/33/05/029 Bajour (7/11/2005) 

— NW/33/05/014 Bajour (2/8/2005) 
NW/30/05/001 Peshawar (11/1/2005) 
I — NW/61/06/036 Bannu (10/7/2006) 
L NW/61/06/019 Bannu (25/5/2006) 

NW/64/06/006 N.Waziristan (9/3/2006) 

Jl— NW/61/06/016 Bannu C2 (18/5/2006) 

NW/62/06/024 LakkiMarwat (26/8/2006) 

NW/61/06/024 Bannu (28/7/2006) 

NW/61/06/027 Bannu (25/7/2006) 

I NW/64/06/008 N.Waziristan (13/9/2006) 

I NW/34/06/021 Khyber (20/10/2006) 

I AFG/06/05/375 (846)C4 
^AFG/06/05/375 (845)C3 

NW/34/07/036 Khyber (22/12/2007) 

I NW/40/07/086 Mardan (16/12/2007) 
~l NW/40/07/082 Mardan (12/11/2007) 
NW/40/07/089 Mardan (21/12/2007) 
I — NW/24/07/800 (Australia) 27/6/2007 
NW/40/07/088 Mardan (15/12/2007) 



Lineage h 



Lineage i 



0.005 

Figure 3 Dendogram of Cluster B summarizing three lineages (g-i) sequence relationship among wild typel polioviruses isolated in 
Pakistan, 2005-2007 based on VPl nucleotide sequences. Tree was constructed by Neighbor-Joining method using Kimura2-parameter 
model. Isolates are indicated by epid no, date of onset and province codes. Codes were mentioned as PAK: Pakistan, AFG: Afghanistan, BN: 
Balochictan, KP: Khyber Pakhtoonkahw, PB: Punjab, SD: Sindh. 
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Laghman reflecting the cessation of this lineage in Paki- 
stan but sharing the same infection with Afghanistan. 

The genetic relationship among the viruses categor- 
ized under lineage 'i' (Figure 3) revealed its prevalence 
in four districts of Pakistan (Khyber, North Waziristan, 
Bannu and Lakki Marwat) as well as Nangarhar district 
of Afghanistan during 2005. But during 2006, it was 
restricted in district Mardan only. In 2007, no infection 
caused by viruses from this lineage except only one wild 
typel poliovirus isolate; a case of international public 
health concern (NW/24/07/800) of importation to Aus- 
tralia by 22 years adult who acquired poliovirus infec- 
tion during his visit to Swat district [38]. 
Cluster C 

This cluster comprised only one isolate from district 
Sanghar (Sindh) which has sequence similarity with the 
2004 wild type 1 polioviruses (Figure 5). The correlation 
proximity of SD/32/05/001 isolate revealed a continuous 
unchecked low level of circulation in this district for 
many previous years. However, no subsequent viruses 
have been reported from this group during 2006-2007 
(Figure 1). 
Orphan lineages 

The covert circulations of 11 wild orphan polioviruses 
(isolates with > 2% sequence diversity with their closely 
related isolate) have been observed in several districts/ 
provinces of both countries (Table 4). During 2005, 
seven orphan lineages were detected mainly from all 
four provinces of Pakistan. No such orphan virus has 
been detected during the same period from AFG. How- 
ever, in 2006 two orphan lineages detected in D.I. Khan 
district of Pakistan and Kandhar province from Afghani- 
stan were grouped in Cluster A. Only two orphan 



- PAK04-2005722999 PAK/PB/22/04/016 Hafizabad 
PAK04-2004737844 PAK/PB/71/04/016 Bahawalpur 

— PAK04-2005720946 PAK/SD/77/04/007 Karachi Saddar 
PAK04-2005717467 PAK/SD/2 1/04/044 Hyderabad 

-PAK05-033 PAK/SD/32/05/001 Sanghar 
-PAK04-2005717767 PAK/SD/31/04/033 Mrpurkhas 
-PAK04-2005717464 PAK/SD/31/04/029-C5 Mrpurkhas 
-PAK04-2004703830 PAK/SD/3 1/04/020 Mrpurkhas 
-PAK04-2005717460 PAK/SD/3 1/04/029 Mrpurkhas 

-PAK04-2005717462 PAK/SD/31/04/029-C4 Mrpurkhas 



Figure 5 Dendogram of Cluster C summarizing sequence 
relationship among wild typel polioviruses isolated in 
Pakistan 2004-2007 based on VP1 nucleotide sequences. Tree 
was constructed by Neighbor-Joining metliod using Kimura2- 
parameter model. Isolates are indicated by epid no and province 
codes. Codes were mentioned as PAK: Pakistan, AFG: Afghanistan, 
BN: Balochictan, KP: Khyber Pakhtoonkahw, PB: Punjab, SD: Sindh. 



viruses have been detected during 2007, one from thatta 
district of Sindh and second from Khyber agency in KP. 

Discussion 

The data presented here illustrates molecular characteri- 
zation of wild type 1 polioviruses endemic in Pakistan 
and Afghanistan during last 3 years, 2005-2007. The 
sequence of VPl gene was employed to determine the 
epidemiological links among the wild type 1 polioviruses 
circulating across the both countries. The VPl gene com- 
prised of 906 nucleotide bases with sufficient heteroge- 
neous motifs that clearly discriminate between various 
genotypes. In 1993 Huovilianen et aL, has been detected 
multiple genotypes of wild type 1 polioviruses in Karachi 
district of Sindh province in Pakistan [39]. However, in 




Figure 4 Geographical location of wild typel polioviruses 
isolated during 2005-2007 from Pakistan-Afghanistan. 



Table 4 Orphan Poliovirus isolates with equal to or more 
than 2% nucleotide sequence divergence from their 
parent isolate reported during 2005-2007 in Pakistan 



Isolate name 


District/province 


Year 


% Divergence 


NW/63/05/006 (2608) 


Tank/KP 


2005 


2.1% 


NW/30/05/001 (081) 


Peshawar/KP 


2005 


3.2% 


NW/30/05/098 (5373) 


Peshawar/KP 


2005 


2.1% 


BN/1 4/05/01 5 (5805) 


K.Abdulah/Balochistan 


2005 


2.3% 


BN/24/05/003 (3845) 


Musakhel/Balochistan 


2005 


2.8% 


SD/5 1/05/037 (4937) 


Larkana/Sindh 


2005 


2.2% 


PB/61/05/005 (225) 


Multan/Punjab 


2005 


2.8% 


NW/60/06/002 (743) 


D.I.Khan/KP 


2006 


2.3% 


AFG/08/06/463 (125) 


Kandhar 


2006 


2.8% 


SD/22/07/019 (4026) 


Thatta/Sindh 


2007 


2.8% 


NW/34/07/036 (6212) 


Khyber/KP 


2007 


2.4% 



The isolate name indicates (Province code/District code/Year of paralysis 
onset/Serial No. of case from the same district) 
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Pakistan and Afghanistan only one genotype (SOAS, 
South Asian) has been endemic for many previous years 
in contrast to African countries where endemicity of 
multiple genotypes have been reported [40,41]. 

Despite this monophyletic character, distinctive reser- 
voirs of wild type 1 polioviruses have been observed that 
confined to specific geographic regions in both coun- 
tries. Based on sequence diversity of < 5% among the 
VPl nucleotide sequence, the wild type 1 polioviruses 
have been classified into certain groups called 'clusters'. 
A phylogenetic tree was reconstructed on the genetic 
relationship of VPl gene sequence for the wild polio- 
viruses detected during the 3 year period which clearly 
displayed the distribution of isolates among three clus- 
ters (A, B, C) (Figure 2). Within each clusters a definite 
pattern of virus distribution was found, representing a 
different lineage or transmission route in different geo- 
graphic regions of Pakistan and Afghanistan that also 
highlights the epidemiological links during their evolu- 
tion from 2005 to 007 time intervals. 

The molecular data obtained from this study explicitly 
reveals the transmission pattern of wild type 1 polio- 
viruses circulation in Pakistan and Afghanistan as a suc- 
cessful tool to monitor the AFP surveillance as well as to 
formulate the mop-up immunization activities estab- 
lished as a reference standard policy to interrupt the wild 
poliovirus multiplicity within population. The underlying 
factors posing a incessant threat to achieve the polio-free 
country status are the areas with poor security and hard 
to reach children communities leaving a number of chil- 
dren unimmunized despite a number of supplementary 
immunization campaigns in the country. Furthermore, 
parallel circulation of some lineages/multiple chains 
among both countries highlights the ongoing continuous 
presence of wild polioviruses. 

The factors behind this worrisome situation are attrib- 
uted to the security concerns and also to the unchecked 
cross-border movement as well. The similar concern has 
already been shown by Kew., et al. [26]. However, the 
synchronized immunization activities across the border 
will be helpful to stop the importation of wild polio- 
viruses from each side of the border. Because immuniza- 
tion has been proved to be a successful approach for 
elimination of poliomyelitis by locating and maintaining 
a high level of immunity in children less than 5 years of 
age. Therefore extended use of type 1 monovalent oral 
poliovirus vaccine (mOPVl) and type 3 monovalent oral 
poliovirus vaccine (mOPV3) selected in areas of Pakistan 
and Afghanistan was adopt to knockdown the virus in 
2005 and 2007 respectively [42-44]. But transmission of 
virus could not be stopped in Khyber, South Waziristan, 
and the areas due to war and armed conflict which 
destroyed almost all public health infrastructures and 



health-care workers could not safely deliver vaccine and 
active surveillance was very difficult to carry out. 

Current molecular data also explores the evolutionary 
pattern of wild type 1 poliovirus isolates in relation to 
their ancestral descent. Based on the 1% (in vivo) nucleo- 
tide evolutionary rate of VPl gene, any wild poliovirus 
isolate indicating equal to or more than 2% nucleotide 
sequence divergence from its ancestor is recognised as an 
"orphan virus" [45]. Findings of this study also put 
emphasis on some silent multiple independent chains of 
transmission (orphan lineages) over a wide geographical 
areas involving cross-border population movement 
between Pakistan and Afghanistan. Regardless of the 
view that orphan viruses reflect some collapse of surveil- 
lance even with superfluous efforts reporting or long 
term circulation without causing overt cases of polio and 
gaps in immunization, the efforts of health care person- 
nel working towards the Polio Eradication Initiative (PEI) 
under the prevalent difficult situation must be 
acknowledged. 

Moreover, epidemiological records and results of this 
study clearly highlights that transmission of wild type 1 
polioviruses is still uncontrolled. However, significant 
progress has been made to localize the wild poliovirus. 
The number of infected districts has been significantly 
decreased from 2005 to 2007 (Table 2) in both countries 
despite of war conflicts, poor access and a number of 
other reasons. In order to achieve polio free status, 
interruption of wild poliovirus transmission is the key 
target. Pakistan is under strict apprehension of policy 
implementation constrains required for high level popu- 
lation immunity. A concomitant locally appropriate 
communications and social mobilization strategy must 
be pursued to improve public awareness for effective 
oral polio vaccine coverage to susceptible population to 
ensure substantial progress. 

Conclusion 

In conclusion, this study provides requisite support to 
find out the locations, the extent of wild typel polio- 
virus circulation in endemic areas, identify reservoir 
communities sustaining wild poliovirus endemicity, 
genetic relationships among these isolates and the 
source of imported wild poliovirus in Pakistan and 
Afghanistan. Moreover these analyses are valuable for 
monitoring the AFP surveillance and to target supple- 
mentary immunization activities with oral polio vaccine 
in order to interrupt chains of virus transmission. 
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